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We present a pliase field crystal model for driven systems which describes competing effects 
between thermally activated difTusional processes and those driven by externally imposed ballistic 
events. The model demonstrates how the mesoscopic Enrique and Bellon [Phys. Rev. Lett. 84, 
2885 (2000)] model of externally induced ballistic mixing can be incorporated into the atomistic 
phase field crystal formalism. The combination of the two approaches results in a model capable of 
describing the microstructural and compositional evolution of a driven system, while incorporating 
elasto-plastic effects. The model is applied to the study of grain growth in nanocrystalline materials 
subjected to an external driving. 



I. INTRODUCTION 

There has been an explosion over the past couple of 
decades in the investigation of nanocrystalline (NC) ma- 
terials as viable replacements for components in vary- 
ing industries, from microelectronics and integrated cir- 
cuits to core components in nuclear reactors. Due to the 
small size features exhibited by NC materials, and their 
propensity to have a high density of grain boundaries 
and phase boundaries; the mechanical, chemical, electri- 
cal and optical properties of NC materials differ quite 
markedly from their polycrystalline counterparts. 

Whether for tuning nanostructure size features in com- 
ponents during their fabrication (microelectronics) or as 
full size components in nuclear reactors, the microstruc- 
tural and phase stability of these components, when ex- 
posed to external driving, is paramount. These external 
processes may be attributed to, neutron or electron irra- 
diation (in nuclear systems) , ion-beam mixing (in fabrica- 
tion of microelectronics) or continued plastic deformation 
(leading to cyclic fatigue, e.g. high energy milling) Of 
importance in these systems is the structural integrity of 
components. The system is driven out of equilibrium (as 
described by thermally activated processes, i.e., thermo- 
dynamic equilibrium), and in the thermodynamic limit, a 
system may exhibit chemical or structural modifications 
manifested in microstructural and phase alterations. The 
resulting long-time configurations of these systems are 
then dependent on the conditions responsible for exter- 
nal driving. 

For polycrystalline materials, there are several exhaus- 
tive reviews and no shortage of literature on materials 
under external driving conditions, particularly when the 
forcing condition is a flux of energetic particles (such 
as electrons, ions or neutrons) Several interest- 

ing and unique phenomena can occur in these materi- 
als. Examples include precipitation in under-saturated 
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solid solutions [4], the lowering of the critical point in 
second order order-disorder transitions d i], the pat- 
terning of dislocations |6| and the amorphization of crys- 
talline materials 0] . Note that external driving is also 
usually accompanied by the production of excess point 
defects, which among other things enhance thermal dif- 
fusion. The common denominator in all driven systems, 
is that they achieve an altered dynamical state, where the 
configuration of such a state is a function of the mate- 
rial properties and forcing conditions The resulting 
state can be rationalized as being caused by the long- 
time effect of two parallel yet competing mechanisms: 
thermally activated atomic jumps and the forced (ather- 
mal ballistic) atomic jumps arising from external forcing. 
Research suggests that a method to alleviate the effects 
of these athermal ballistic effects is the use of NC ma- 
terials, which feature a very high density of potential 
point defect sinks such as grain boundaries or interphase 
boundaries. Recently, Bai et al. fsl have outlined a theo- 
retical explanation for the apparent resistance of NC ma- 
terials (pure Cu) to ballistic conditions. This resistance 
to ballistic effects has also been observed experimentally 
by Rose and coworkers 01 on NC Pd and Zr02 mate- 
rials, Hochbauer et al. |lO| on Cu-Nb multilayers and 
in various atomistic simulations [111, [12] • Unfortunately, 
most existing atomistic models are not capable of cap- 
turing the long-time scales relevant in studying sustained 
ballistic effects. While techniques have been developed 
to accelerate atomistic simulations Q, it is not clear if 
these methods can be extended to a time scale appropri- 
ate for the solid state diffusion mechanisms controlling 
the microstructural and phase stability of NC materials. 

Enrique and Bellon [l^ have recently developed a 
continuum description of alloys subjected to the exter- 
nal driving via ballistic events. The underlying theme 
of their model is the competing dynamics between two 
dynamical mechanisms: one which is thermally driven to 
bring the system to thermodynamic equilibrium; and the 
other is athermal particle exchanges driving the system 
out of equilibrium, an idea which dates back to Mar- 
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tin [l5|. This is accomplished through an effective free 
energy functional of the Cahn-Hilliard form, whereby 
they forgo the explicit description of discrete interac- 
tions and supersaturation of excess point defects and in- 
stead make use of effective interactions to account for 
the ballistic contributions. The model has been used to 
study compositional patterning in alloys driven by irra- 
diation [14i, illl, phase stability of alloys under irradia- 
tion [l3| , non-equilibrium fluctuations in alloys under ir- 
radiation [l3| and irradiation induced spinodal decompo- 
sition in the presence of dislocations [18:] . The model can 
capture the appropriate long-time scales needed in ex- 
amining sustained ballistic effects, unfortunately it does 
not incorporate different crystalline orientations, grain 
boundaries, defects or elasto-plastic effects. These are 
important in the microstructural and phase stability of 
NC materials under sustained external forcing. 



In recent years a new atomistic modelling paradigm 
coined the phase field crystal (PFC) method has 
emerged. This approach is an effective modelling tool 
that describes microstructure evolution over diffusive 
time scales, while incorporating atomistic-scale elastic- 
ity, topological defects and grain boundaries. Several 
PFC models have been developed and used successfully in 
the description of solidification [20j , spinodal decomposi- 
tion [2l|, elasto-plasticity [l^], thin film growth and 
structural phase transformations in pure materials (2^ 
and alloys (25j . The purpose of the present work is to 
incorporate the mesoscale model of ballistic effects in 
driven systems developed by Enrique and Bcllon fli*] into 
the PFC formalism to create a simulation framework ca- 
pable of describing the evolution of composition and mi- 
crostructure under ballistic conditions. Like the original 
model, the present model does not explicitly treat the 
supersaturation of point defects, and instead treats the 
long-time ensemble average of the discrete ballistic effects 
as effective interactions. Our aim is to provide a mod- 
elling framework having the capabilities of simulating the 
long-time scales of sustained external forcing, while in- 
corporating the atomic effects mentioned previously. We 
will demonstrate the properties model by examining as 
a first example the grain growth behaviour in the PFC 
formalism of NC materials under forced ballistic condi- 
tions. 



The paper is organized as follows. First, we present the 
PFC model and its modification to account for ballistic 
effects. After introducing the model, we analyse the long 
wavelength properties of the model, where we highlight 
modifications to the model of Enrique and Bellon. We 
follow with numerical simulations of grain growth under 
ballistic effects. Finally we conclude and summarize. 



II. EFFECTIVE PFC ENERGY FUNCTIONAL 
FOR BALLISTIC MIXING 

We begin with the standard binary alloy PFC energy 
functional in scaled form [21| , 
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where the energy is scaled by kBTpiR'^. The field n is 
the dimensionless local density, ijj the dimensionless con- 
centration, kg is the Boltzmann constant, T the temper- 
ature, Pi the reference liquid density, and R the average 
atomic radius, with d being the dimensionality. Note 
the final three terms in Eq. ([1]) comprise a Cahn-Hilliard 
equation [il], which has been used to study phase sep- 
aration. B^ and i?^ are the dimensionless bulk moduli 
of the liquid and solid, respectively, setting the energy 
scale of the system. Following Ref. [13], here we con- 
sider the expansion of the moduli where Bi = Bq + B2V'^ 
and B^ = Bq. The remaining parameters are constants, 
which in principle are functions of the direct two point 
correlation functions, and can be calculated from first 
principles or fit to phenomenological databases or theo- 
ries of materials properties, e.g. surface energy [27l - [29| . 
However, the current work will only treat these param- 
eters as constants. In the liquid phase n is constant ev- 
erywhere, while in the solid it assumes a spatially peri- 
odic structure. In two dimensions (2D), the functional in 
Eq. ID) has a phase diagram of coexisting liquid and solid 
phases with hexagonal symmetry. For greater detail in 
the derivation of the PFC energy function presented in 
Eq. H]), the reader is referred to Ref. [21*]. 

The density and concentration fields, n and ip, follow 
dissipative dynamics that minimize J-. To include exter- 
nal forced ballistic effects, the dynamics are augmented 
by two Enrique and Bellon like source terms in the den- 
sity and concentration equations. They are written as 
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with M the mobility and F the frequency of forced atomic 
exchanges, proportional to the flux of incident particles. 
Note that Eqs. ^ and ([3]) are deterministic, however 
to incorporate the short time fluctuations of the system, 
in Section llVI when we report on numerical simulations, 
we append to these equations two stochastic variables, 
C and ^, having a Gaussian distribution with zero mean 
and amplitudes Aq and for the density and concen- 
tration flelds, respectively. {ip)-fi and (n)7j denote the 
corresponding weighted spatial averages of the concentra- 
tion and density flelds respectively, in response to forced 
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atomic exchanges with an average distribution distance 
TZ. The weighted averages are defined by, 



-r')^(r') 



and 



{n)n^ / (ir'wK(r-r')n(r'). 



(4) 



(5) 



The function w-n is a weight function describing the long- 
time spatial extent of ballistic exchanges in crystalline 
materials. A Yukawa potential was chosen as a form for 
the weight function in Ref. , for the case of exchanges 
driven by irradiation. However in the present work we 
adopt the 2D analog; 
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where Kq is a modified Bessel function of the second 
kind. As stated above, TZ describes the average spatial 
extent of ballistic events and is also related to the energy 
of the ballistic driving force. We note that the above is 
the lowest order description of ballistic events in driven 
systems that can be included for a binary alloy model 
that has coupling between density and composition vari- 
ations. In general T can take on different values for the 
density and concentration field, respectively, and two sep- 
arate functions and can be introduced for each re- 
spective field. The model presented here, like the model 
of Enrique and Bellon, has the underlying idea that a 
system exposed to external driving can be described by 
two mechanisms acting in parallel: thermal diffusion de- 
scribed by terms of the formMV^T^ and ballistic events 

on 

described by terms of the form F (n— {n}Tz). Similarly 
for the concentration order parameter. Note that when 
TZ ~ 0, dynamics are driven only by the thermal portions 
of Eqs. @ and ©. 

Including only the lowest level of description, we con- 
struct an effective PFC energy functional for ballistic 
events, which affords an opportunity to explore equilib- 
rium properties under ballistic conditions. Following En- 
rique and Bellon, we can write the effective functional 
as E ^ T + 7®, where (5 is the external energy source 
representing effective interactions arising due to ballistic 
effects and 7 = F/M, the reduced forcing frequency. The 
effective interactions are described by a self interaction 
energy written as 

© = iy"drdr' {n(r)g(r-r')n(r') + ip{r)g{r-r')ip{r')} , 

(7) 

where the kernel for ballistic exchanges g satisfies the 
equation V'^g{r — r') — — [S{r~r') — WTz{r — r')]. Which is 
to say, the kernel satisfies the Poisson solution to a point 
source perturbed by a weight function. 



III. LONG WAVELENGTH AND EQUILIBRIUM 
PROPERTIES 



We first interrogate the long wavelength (i.e., phase 
field limit) limit of the model. Secondly, we examine 
its bulk properties under forced driving, i.e., a kinetic 
phase diagram under ballistic conditions. To examine 
the long wavelength properties of the model, we perform 
an amplitude expansion of our effective PFC free energy, 
E. To perform the expansion, we will be substituting 
into the effective free energy a single mode expansion for 
the density. 



?^o(r) 



-|- c.c. 



(8) 



7io(r) is the dimensionless average density of the sys- 
tem, a conserved quantity which is slowly varying spa- 
tially with respect to Gj, r]j{r) are dimensionless spa- 
tially varying complex amplitudes which like the average 
density are slowing varying on atomic scales, Gj repre- 
sents the lowest set of reciprocal basis vectors necessary 
to describe the crystal structure of interest and c.c. de- 
notes the complex conjugate. The result is then coarse 
grained using the volume averaging method of Refs. [28j 
and [2^. For a crystal of hexagonal symmetry in 2D, 
coarse graining E gives, to lowest order in the ballistic 
terms; 

^ = / dr i^J2{ABo+r9i\Gjl'R-) + Siynl-2tno) 
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where A£?o = Bq — Bq defines a temperature scale, the 
caret denotes the Fourier transform, with g = TZ^/{i + 
|GJ|27^2) and Gj = -f 2iGj ■ V. The free energy func- 
tional of Eq. ^ describes a system with multiple crys- 
tal orientations, elasto-plastic effects through the com- 
plex nature of the amplitudes, as well as density changes 
through the field rio- Moreover, as a coarse grained 
model, it operates on scales much larger than the lat- 
tice constant of the solid. A noteworthy result of Eq. ^ 
is that the "fg{\Gj\,TZ) term competes with the temper- 
ature (A_Bo), as both are modulated by the amplitude 
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magnitude squared, \r]j\'^. In other words forced external 
driving is manifested, partly, as shifts in the tempera- 
ture scale leading to a new effective temperature. This 
is in agreement with previous theories of phase equilib- 
ria under driven conditions 0, IT5| . Furthermore, the 
ballistic term in concentration, 7^ / dr'g{r — r')'ip{r'), 
will be shown to renormalize the coefficient of the ip'^ in 
the Cahn-Hilliard portion of the energy in the thermody- 
namic limit, i.e., limit TZ — oo. Namely, the model also 
predicts that ballistic effects cause alterations in the crit- 
ical transition temperature as well. We also note that in 
the absence of density variations, defects and the elasto- 
plastic effects inherent in the complex nature of the am- 
plitude formulation, i.e., all |?7j| = const, and scale out of 
the problem, we recover the Cahn-Hilliard ballistic model 
of Enrique and Bellon [IJ]. The repercussions of these 
effects are explored below when we construct phase dia- 
grams for this model. 

A. Pure Material 

The equilibrium phase diagrams follow standard mini- 
mization procedures analogous to Refs. and ^2T\. For 
a pure material, we assume a real and constant amplitude 
rjj = (p for all j and a constant wave vector, to describe 
the periodicity of the solid, ignore all terms in -0 and set 
t — and v — 1. Making a long wavelength approxi- 
mation, i.e., qTZ 3> 1, we first begin by minimizing the 
energy for the equilibrium wave vector, q^q — \Gj\, which 
is substituted back into the energy. Next, the energy is 
minimized for the real amplitude (f>. After substitution of 
the minimized amplitude back into the energy, we have 
a resulting energy which is a function of the conserved 
quantity rio, temperature (ABq), and the reduced forc- 
ing frequency 7. A common tangent or the equivalent 
Maxwell equal area construction leads to Fig. [TJ which 
plots the a dynamical phase diagram in ABq, rio and 7 
space. We notice a trend of decreasing solid-liquid tran- 
sition temperature with increasing 7. As previously dis- 
cussed, we can attribute this behaviour to the competi- 
tion between the temperature and ballistic term, where in 
the thermodynamic, constant amplitude limit, the com- 
petition term takes the form of (Ai3o + l/q)<t)'^ ■ 

B. Eutectic Binary Alloy 

Moving on to the binary alloy, we again assume real 
and constant amplitude rjj = cf) for all j and a constant 
wave vector, q, to describe the periodicity of the solid. 
We set = -1.8, Bg^ = 1, t = 0.6, = 1, m = 4 and 
rio = 0- A spinodal and a eutectic alloy are differen- 
tiated, in this PFC model, by choice of the parameter 
(jj. Here we report the calculation for the eutectic alloy 
only, i.e. lu = 0.008. After the parameters are set, we 
follow the same minimization steps as outlined for the 
pure material in the preceding section with respect to 
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FIG. 1. (colour online) Kinetic phase diagram for a pure ma- 
terial exhibiting density variations in ABq x Uo x "f space, 
shown here projected onto the ABo x Uo space for choice 
values of the reduced forcing frequency 7. Squares with 
lines(blue) represent the standard PFC model, 7 = 0, dashed 
lines (green) correspond to a ballistic case of 7 = 0.02, while 
squares (black) correspond to a ballistic case of 7 = 0.05. 



and (jcq. With the average density being chosen to be 
zero, the resulting minimized energy is a only a function 
of '0, temperature (Ai?o), reduced forcing frequency 7 
and ballistic distance TZ. A common tangent construc- 
tion or equivalent Maxwell equal area construction gives 
Fig. In Fig. [21 we present the phase diagram of a eu- 
tectic alloy under ballistic effects in ABq, V' ^^id 7 phase 
space. In the figure, a = 27r/geq is the lattice spacing 
associated with the hexagonal unit cell in 2D. There are 
several things to note. There is still a general shift of 
the solid-liquid transition temperature, additionally how- 
ever, we also witness the alteration of the eutectic alloy 
into a spinodal alloy. Again, in the thermodynamic limit 
the coefficient of the amplitude square term becomes 
(ABo + 7/9)0^7 which from the pure material calcula- 
tion above, was shown to affect the solid-liquid transition. 
This term is what primarily leads to the general shift in 
the solid-liquid transition. However, for the binary alloy, 
the coefficient multiplying the ip"^ in the Cahn-Hilliard 
contribution becomes (w + 125^0^ -I- ^'R?) ^p^. The jTZ'^ 
(i.e. qTZ S> 1) contribution to this term simultaneously 
accounts for the change from eutectic to spinodal, the 
ever increasing range of the solid solution region and con- 
tributes to the change in the solid-liquid transition lines 
in our model as well. 



IV. DYNAMICS 

Two sets of numerical simulations were conducted to 
study the role of ballistic driving on the grain growth 
of a pure nanocrystalline material. The first is designed 
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7 = -- 7 = 0-02 K = 2a • • 7 = 0.04 K = 2a 




FIG. 2. (colour online) Binary eutectic phase diagram under 
forced ballistic driving, circles and lines(black) represent the 
standard binary PFC model, 7 = 0, dashed lines (red) cor- 
respond to a ballistic case of 7 = 0.02 and TZ — 2a, while 
circles(blue) correspond to a ballistic case of 7 = 0.04 and 
TZ — 2a with a being the lattice spacing. 



the fastest coarsening rates. 

The results of Fig. [3] can be understood as follows. 
The long-time ballistic effect acts over a length scale TZ 
and within this spatial extent the mobility and/or driv- 
ing force is enhanced. Thus, when TZ becomes compara- 
ble to the size of grain "4" , the coarsening rate is more 
pronounced as seen in the TZ — 5a and 10a cases. The 
plateaus evident at late time for these two conditions are 
a result of the background signal of the spectrum after 
the disappearance of grain "4" from the system. When 
the ballistic distance is much larger than the size of grain 
"4" , i.e., TZ = 50a and 200a in the second regime, there is 
equal enhancement in the competition from all grains for 
growth (although still enhanced compared to the 7 = 
case). This results in a deadlock, which does not allow for 
preferential growth. Concerning the TZ — la, the spatial 
extent of this ballistic event is not large enough to im- 
pact any relevant length scales necessary for pronounced 
preferential growth. However having influence on all in- 
teratomic length scales in the system, particularly those 
associated with boundary widths, it effectively behaves 
like the larger ballistic distances. 



to illustrate the physics of the proposed ballistic PFC 
term. The second set of simulations examines the grain 
growth behaviour of a polycrystalline sample under bal- 
listic conditions. For the simulations, we solve Eq. 
with Ax = 0.785, At = 1, Af = 1, ABq = -0.26, 
Uo = 0.285 and C 7^ (with amplitude — 0.01), using 
a semi-implicit Fourier technique. 



A. Coarsening of a 3-Sided Grain 

For the first set of simulations, we have 4 grains as 
shown in the inset of Fig.[3]as an initial structure. In this 
configuration, conventional theory states that grains "1" , 
"2" and "3" will tend to grow at the expense of grain "4" , 
driven by gradients in the chemical potential resulting 
from curvature effects. Simulations were performed for 
the regular PFC dynamics (7 = 0), for the proposed 
ballistic PFC model at 7 = 0.01 and for several values 
of the ballistic distance, i.e. TZ — la, 5a, 10a, 50a and 
200a. For a particular measure of the coarsening rate, 
we focus on the behaviour of grain "4". Results were 
averaged over several runs, where grains "1" , "2" and "3" 
assumed different orientations while keeping grain "4" 
fixed. In Fig. [31 we plot the average intensity of the Bragg 
peaks obtained from the power spectrum for grain "4" 
(which is proportional to its area) versus time. Several 
interesting aspects of the plot are noteworthy. First, the 
rate of coarsening is enhanced for all cases of ballistic 
mixing. However there exists two regimes of noteworthy 
behaviour. Namely, we find one regime wherein the TZ = 
5a and 10a conditions cluster and another for the TZ = 
la, 50a and 200a conditions. In the first regime we have 
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FIG. 3. (colour online) Grain intensity versus time for a 
512Aa; x 512Aa;, with 7 = 0.01 for varying varying TZ. Left 
bottom inset shows an example of the initial microstructure 
used to perform the simulations. Density (n) map is shown, 
where lighter regions represent areas with a higher probabil- 
ity. 



B. Grain Growth 

The second set of simulations examines grain growth 
for a nanocrystalline sample in a 1024Ax x 1024Ax sys- 
tem, with approximately 100 grains. Figure |3] plots the 
average grain size versus time, for 7 = 0.05 for vary- 
ing ballistic distances TZ. Grain growth was quantified 
by Fourier analysis, where the full width at half max- 
imum of the the Bragg peak was measured in the ra- 
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dially averaged signal. The multigrain simulations are 
consistent with the results of Fig. [31 Grain growth is en- 
hanced for all TZ values. Also note, in the time range of 
500 — 1000, the average grain size is 75a, and the growth 
is enhanced most strongly for TZ comparable to this grain 
size [TZ = 100a). Therefore, our results demonstrate 
that nanocrystalline materials, although possibly offer- 
ing increased resistance when exposed to ballistic driving 
forces, such as that from irradiation damage, are suscep- 
tible to enhanced grain growth under such driving forces, 
with the rate of growth generally increasing with increas- 
ing ballistic energy (i.e. TV), and particularly enhanced 
when TZ is comparable to the grain size. 




time 



FIG. 4. (colour online) Grain growth results plotting average 
grain size versus time, on a log-log scale, for a 1024Ax x 
1024Ax. Here 7 = 0.05 for 7^ = 10a, 100a, 200a and 7^ = 
1000a. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have developed an atomistic contin- 
uum model of competing effects between thermally acti- 
vated diffusional processes and those driven by externally 
imposed ballistic events in both pure materials and alloys 
using a PFC free energy functional. As an extension of 
the model of Enrique and Bellon, our approach is capable 
of describing defect microstructure, multiple crystal ori- 
entations, and elasto-plastic effects. A long wavelength 
analysis of our model predicts that ballistic events change 
the phase diagram through a shift in the solid-liquid tran- 
sition temperatures for a pure material. For a eutectic 
alloy, the model further predicts an alteration of equi- 
librium through shifts in the critical transition tempera- 
tures in the solid solution regions, an effect which have 
long been conjectured in the irradiation damage litera- 
ture. Through an investigation of the kinetic proper- 
ties of the model, we find that under ballistic driving, 
grain growth is enhanced with increasing average ballis- 
tic distance, with a predicted enhancement of grain sizes 
comparable to the ballistic range TZ. Further numerical 
and experimental investigation is required to elucidate 
the possible myriad of behavioural aspects in nanocrys- 
talline system when subjected to ballistic driving. 

We acknowledge the National Science and Engineering 
Research Council of Canada (NSERC) for financial sup- 
port and David Monticl and Jonathan Stolle for useful 
discussions. 
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